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Abstract 

Pair creation by a photon in an oriented crystal can be considered in a 
frame of the quasiclassical operator method, which includes processes with 
polarized particles. Under some quite generic assumptions the general ex- 
pression is derived for the probability of pair creation of longitudinally 
polarized electron (positron) by circularly polarized photon in oriented 
crystal. In particular cases t?q -C Vo/m and $0 S> Vo/m (i?o is the an- 
gle of incidence, angle between the momentum of initial photon and axis 
(plane) of crystal, Vo is the scale of a potential of axis or a plane relative 
to which the angle $0 is defined) one has constant field approximation and 
the coherent pair production theory correspondingly. Side by side with 
coherent process the probability of incoherent pair creation is calculated, 
which differs essentially from amorphous one. At high energy the pair cre- 
ation in oriented crystals is strongly enhanced comparing with amorphous 
medium. In appendixes the integral polarization of positron is found in 
external field and for coherent and incoherent mechanisms. 
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1 Introduction 



The study of processes with participation of polarized electrons and photons 
permits to obtain the important physical information. Because of this reason 
the experiments with use of polarized particles are performed and are planning 
in many laboratories (CERN, Jefferson Nat Accl Fac, SLAC, BINP, etc). In this 
paper it is shown that oriented crystal is an unique tool for work with polarized 
electrons and photons 

The quasiclassical operator method developed by authors is adequate 

for consideration of the electromagnetic processes at high energy. The probability 
of polarized pair creation by a circularly polarized photon has a form (see [3], 
p. 73, Eq.(3.12); the method is given also in j^,©) 
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where for polarized electrons and positrons one has within the relativistic accu- 
racy 

fee' 



= -^=<(p)eW(p') (1.2) 



here k v = (o>,k) is the 4-momentum of the initial photon, p M = (e, p) is the 
4-momentum of the created electron, k 2 = 0, x^(t) = (t, r(t)), t is the time, and 
r(t) is the created electron location on a classical trajectory, kx(t) = ut — kr(t), 
e' = u) — e, p' = k — p, A = ±1 is the circular polarization of photon, s(s) = ±1 
is the longitudinal(with respect to direction of motion) polarization of electron 
(positron), we employ units such that h — c — 1. 

Passing on to two-component spinors (p and choosing as the quantization axis 
(the axis z) the direction of photon momentum k we get 



m s + m I Lfs 



Here s± = (a x ± icr y )/2 is the raising (lowering) operator: 

a z ip s = sip s , s x ps = Sx-sf-s, (fith = Sx^-s- 

Below we will consider the small angle approximation (the vector p is directed 
nearly along the vector k, just this configuration yields the main contribution into 
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the process probability) and will retain only the main terms in the decomposition 
over I/7 (7 = e/m is the Lorentz factor). Within this accuracy one obtains 



7 Mp ' ] = T2 
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where i9 = v — n = vj_, v = v(t) is the electron velocity, vj_ is the component 
of electron velocity perpendicular to the vector n = k/u;. Let us note that 
tf' = V - n = p'/e' - n = (n - v)e/e', #7' = $7. 
Using the relations 

ip+a-dips = 0, (p± a tr0(p a = tf x + isfiy = <pfs\(p-s = S\, s 8\,s, 

(p+fTfisxip-s = 6x&pt tr &'Ps = $ {x) 5 x ,s5x-s, (1.6) 
we find within the adopted accuracy 
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Taking advantage of relation 8\ s — 8\ }S = (1 + £)/2, where £ = As, we obtain 
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Here for visualization we have written e = e_ for the electron energy, e' = e + for 
the positron energy, = As, £ + = As for description of electron and positron 
polarization correspondingly, $ 2 ,i = #(£2,1) > 

= (Mh)+iXMtMMti)-i*Mti)) =^i-iX(n(-& 2 x-& 1 )). (1.9) 

Substituting Eqs. (|1.8|) . ()1.9|) into Eq. ljl.lj) one obtains the completely differen- 
tial probability of pair creation for circularly polarized photon and longitudinally 
polarized electron and positron. The above analysis shows that in the situation 
under consideration the derivation with direct calculation of matrix element is 
very simple comparing with the standard procedure (with calculation of traces 
for polarized particles). This is in consequence of Eq. (jl.6|) . 

It follows from Eqs. (jl.7|) . f)1.8j) that in the case when the energy of one of 
particles is essentially larger than other particle energy, then together with energy 
the polarization of photon is transmitted to this particle. This conclusion can be 
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obtained immediately from Eqs. (jl.3|) . (|1.4p if one neglects the angle $ in the wave 
function Eq. (jl.5|) for the particle with large energy and reserves it in the wave 
function for the particle with small energy only Then either one wave function 
or another becomes proportional to either (f s or <^_ 3 . In this case the polarization 
of particle (see Eq. (jl.4j) is defined uniquely. If in this case in addition the angles 
$ ^> I/7 contribute mainly (■$ is the angle between momenta of initial photon 
and created particle), then one can neglect the terms in Eqs. (jl.7J) . ()1.8|) with the 
same polarizations of created particles. Then in Eq. ()1.8p only one term remains 
which defines uniquely the polarizations of both particles. This fact agrees with 
helicity conservation rule in electromagnetic processes. This situation is realized 
in the process of pair creation by a photon in an external field for k>1 jlj, as 
well as for the same process under strong influence of the LPM effect 

Summing in Eq. (jl.8|) over the polarizations of created electron £_ and omit- 
ting the term (v(i?i x $2)) which vanishes at integration over the azimuthal angle 
of created pair we get 



m 



{u 2 (l + £+) + 7 2 2 # 1 [4(1 + £+) + - £+)] } • 

(1.10) 

The cross process of emission of photon with energy uj by an electron with 
high energy e in oriented crystal was considered recently by authors ^2]- It is 
evident that this process possesses the similar properties in the case when the 
photon takes away an essential part of electron energy (e' = e — uj <C e). The 
explicit expression for the transformed combination i2^(i 2 )i2^*(ti) in Eq. ()1.8j) 
can be found using the standard substitutions: 



£_ — > e', e + — > —e, uj — > —u, A — > —A, 



;i.ii) 



As a result we find 
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where £ = AC, £' = AC'- Substituting Eq. fjl.l2|) into the formula for the proba- 
bility of photon emission (see e.g. [I], p. 63, Eq.(2.27)) one obtains the completely 
differential probability of radiation of circularly polarized photon from longitu- 
dinally polarized electron in the case when the final electron has longitudinal 
polarization Summing in Eq. (jl.l2j) over polarizations of final electron C' and 
omitting the term (v(#i x $ 2 )) which vanishes at integration over angles of 
emitted photon we arrive to Eq.(1.8) of 
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It should be noted that while the terms depending on angles in Eqs. (|1.8]) and 
(|1.12|) describe the spin correlations arising from helicity conservation rule, the 
spin correlations in terms oc u 2 describes by rule of conservation of projection 
of angular momentum on the direction of motion at zero angles. This explains 
absence of the term with (1 — £_)(1 — £+) in Eq. (jl.8j) and the term with (1 — 
£)(l + £0mEq.((XI2J). 

It should be noted that a few different spin correlations are known in an exter- 
nal field. But after averaging over directions of crystal field only the longitudinal 
polarization considered here survives. 



2 General approach to pair creation in oriented 
crystal 

The theory of high-energy electron radiation and electron-positron pair creation 
in oriented crystals was developed in [Zj-|H], and given in In these publi- 
cations the process radiation from electron and pair creation by a photon was 
considered for unpolarized particles. Since the expression for the pair creation 
probability has the same structure as Eqs. (|1.8|) . ffl.lOj) below we use system- 
atically the methods of mentioned papers to obtain the characteristics of pair 
creation of longitudinally polarized electron (positron) by a circularly polarized 
photon. 

Let $o be the photon angle of incidence with respect of a chosen axis of 
crystal and V be the scale of the corresponding continuous potential of the axis. 
For i9 S> (Vo/e) 1 / 2 = $ c , created particles are moving high above the potential 
barrier, and can be described in terms of rectilinear trajectory. In rest frame of a 
crystal there is only an electric field of axis E. In frame moving with a relativistic 
velocity v = nv (n = k/a>, k 2 = uj 2 — k 2 = 0) along the photon direction 
momentum, a magnetic field H = 7u(E x v) arises {p/ v = (1 — v 2 )^ 1 ^ 2 ^> 1) and, 
as it is well known, that resultant field in this frame can be represented with 
relativistic accuracy in the form of plane waves with frequencies j v \q\\ \v(q\\ = qn). 

The periodic crystal potential U(r) can be presented as the Fourier series (see 
e.g.jU, Sec.8 ) 

f/(r) = EC(q) e - iqr , (2.1) 

q 

where q = 27r(ni, n^, n$)/l] I is the lattice constant. 

The equivalent photon flux averaged over time and over transverse coordi- 
nates is the sum of partial contributions 3 q . The later quantity has the form 

where q ± = q — n(qn). In the interaction region, the transverse size of which 
is of the order A c = 1/m, and the longitudinal size is the process formation 
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length, which is of the order 2Ti/ r y v q\\ in the center mass frame of the incident 
and equivalent photons, there are N q ~ 2ir\ 2 c \S q /^ v q\\\ photons. The effective 
strength of interaction is characterized by the parameter 



aN, 



ph 



a: 



q 



E 

q 



|g(q)l 

m 2 \q\\ | 



2 



(2.3) 



This parameter is purely classical (it does not contain Planck's constant H) and 
always arises in problems with electromagnetic interaction in external field. For 
otNph <C 1 the external field can be taken into account in perturbation theory, 
while for aN p h ^> 1 one has the constant field limit as it is known from the theory 
of interaction of photon with the plane wave field (see e.g. ||). For estimates one 
can assume |G(q)| ~ V , qn ~ q±$o, in which case 
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ph 
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(2.4) 



Therefore for i9 Vo/m = the constant field approximation is applicable 
while for $o 3> Vo/m the perturbation theory is valid, the first approximation of 
which is the coherent pair production theory (see e.g. [TUj). 

In a crystal one have to integrate over pair creation points. Substituting 
Eq. (|l.lUj) into Eq. (jl.l|) we find for the probability of creation of longitudinally 
polarized positron by a circularly polarized photon 
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where a = e 2 = 1/137, the vector $ is defined in Eq. (|1.5jl . V is the volume 
of the crystal, p + is the momentum of the positron in the creation point ro- 
Corresponding, the probability of creation of longitudinally polarized electron 
by a circularly polarized photon is 
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dt\ / dt 2 e l 
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From Eq. (J2.5|) one can find polarization of created positron 



dvj£ += i — dw. 



dw, 



+ dw, 



*+=-! 



(2.7) 



In the limiting case e + it follows from Eq. (j2.7|) that ( p — > 1. 

The particle velocity can be presented in a form v(t) = v + Av(i), where v 
is the average velocity. If $0 3> $ c , we find Av(t) using the rectilinear trajectory 
approximation for the potential Eq. (j2.1J) 



Av(t) = -- E — Q± exp[-z(?||t + qr)], 



(2.8) 



where q\\ = (qn), q ± = q — n(qn) (for detail see Eqs.(3.27)-(3.31) in jlj). Substi- 
tuting Eq. (j2.8j) into Eq. (|2.5)l and performing integration over u = n — Vo (d 3 p ~ 
£ 2 cfe<iu) and passing to the variables t, r : ii = i — r, £2 = £ + t, we obtain 
after simple calculations the general expression for the probability of creation of 
longitudinally polarized positron by a circularly polarized photon 




r d 3 r 
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where 
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(2.10) 



3 Pair creation for $o Vb/ra 

(constant field limit and corrections to it) 

The behavior of probability Eq. (J2.9|) for various entry angles and energies is deter- 
mined by the dependence on these parameters of the phase A pl given Eq. ()2.10p . 
Here we consider the axial case for t? <C Vq/tti = $ v . The direction of crystal 
axis we take as z— axis of the coordinate system. The order of magnitude of the 
double sum in A v \ is (G/m) 2 (q±/q\\) 2 ^(q\\,q'pT). For the vector q lying in the 
plane (x,y) we introduce notation q t , for such vectors one has q z = and the 
quantities in Eq. (j2.1()j) can be estimated in the following way: 



G(q) 



Vo, 



q± ~ l/a, q\\ ~ -& /a, 



(3.1) 
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where a is the size of the region of action of the continuous potential. For all 
remaining vectors q± ~ q» ~ 1/a. Then the contribution to the sum of the 
terms with q z ^ will be ~ (V /m) 2 ^ < (V /m) 2 . Since (V /m) 2 1 this 
contribution can be neglected. Thus, we keep in the sum only terms with q 4 
for which its value is ~ (Vb/m^o) 2 ^- The large value of the phase A p i leads 
to an exponential suppression of probability dW^ + . Therefore the characteristic 
value of the variable r in the integral Eq. ()2.10|) (which have the meaning of the 
formation time (length) of the process) will be adjusted in a such way that the 
large factor (Vo/m-# ) 2 will be compensated by the function ^/(q^qLr), i.e. for 
small entry angles the contribution gives region where q\\T <C 1. Expanding the 
phase A p i correspondingly we find an approximate expression for $ $ v 
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here f? = rot- We can rewrite Eq. ()3.2j) in the terms of the average potential of 
atomic string U(g) = G(q t ) exp(—iq t g): 
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where b = W(g)/m, V = d/dg. For the pre-exponential factor in Eq. ()2.9|) we 
find 



b 2 + _( b (nV)b) 



(3.4) 



Taking the integral over r we obtain the spectral probability for $o *C Vo/m. 
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where 
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here A = — , K V (X) is the modified Bessel function (McDonald's function). 

3£_£: + |b| 

Since the expression for dW^ + is independent of z, it follows that / d 3 ro/V — > 
/ d 2 g/s, where S is the transverse cross section area per axis. The term in 
Eq. (|3.5|) with Ro(X) represent the spectral probability in the constant field limit. 
The other terms are the correction proportional t9q arising due to nongomogeneity 
of field in crystal. 

If the potential U(g) can be assumed to be axially symmetric, we put U = 
U(q 2 ) and one can integrate over angles of vector g. We obtain 
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where we have gone over to the new variable x = g 2 /a 2 , x < x , Xq 1 = 7ta 2 dn a = 
ira 2 /S, a s is the effective screening radius of the potential of the string, n a is the 
density of atoms in a crystal, d is the average distance between atoms of a chain 
forming the axis. The notation U'(x) = —Vog(x) is used in Eq. ()3.7j) and 



2u 2 V lj dU to /— / \ 

A = --, k s = — ^— , k(x) = — r = 2K s y/xg{x). (3.8) 



For specific calculation we use the following expression for the potential of axis: 
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X +7) I \ Xq + T] 



g{x) 



(X + 7])(x + T] + 1) ' 

(3-9) 

For estimates one can put Vo — Ze 2 /d, rj ~ 2u\/a 2 = 771, where Z is the 
charge of the nucleus, u\ is the amplitude of thermal vibrations, but actually the 
parameters of potential were determined by means of a fitting procedure using 
the potential Eq. ([2.1|) (table of parameters for different crystals is given in Sec. 9 
of 0]). The function k(x) vanishes at x = as a result of thermal vibrations and 
reaches a maximum at 

x = x m = ~ {[1 + 1677(1 + V )] 1/2 -I-277} (3.10) 

and then fall off as x~ 3//2 . It should be noted, that for all crystals 77 -C 1 and in 
this case x m ~ 77 and n(x m ) = K m ~ k s / ^/fj. 

In Fig.l the spectral probability of coherent pair creation dw^ + /dy in tung- 
sten, T=293 K, axis < 111 > is given. The calculation was performed using 
Eq. ([3.7)1 with term Rq(X) only (the corrections oc i?q were omitted). The sum of 
curves at the indicated energy gives unpolarized case. At energies u> = 100 GeV 
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and 00 = 250 GeV the coherent process dominates over the Bethe-Heitler mech- 
anism. 

Let us consider the case of small values of parameter k. In the constant field 
approximation (the term with Rq(X) in Eo. fl3.7j) ) in the probability dW^ + (e + ) 
one can substitute the asymptotic of McDonald's functions K u {\) at A 3> 1. 
After this one obtains using the Laplace method in integration over the transverse 
coordinate x 



ds-L. 
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where K m = K,(x m ), K" m = K,"(x m ). For unpolarized particle this result coincides 
with Eq.(12.14) of 4J. If n m <C 1 one can integrate Eq. (j3.11j) also using the 
Laplace method. The integral probability of polarized pair creation is 
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At relatively low photon energies where the parameter K m is small 

uVq 



(3.12) 



K M 



y/rj yfql V2m 3 m 



V2 



(3.13) 



and the main contribution to e + e~ pair creation gives the incoherent (Bethe- 
Heitler) mechanism. The influence of effective crystalline fields on incoherent 
processes is discussed in jlj, Sections 7, 21, but the polarization effects were not 
included. 

In the case where the influence of effective crystalline field is weak (ki < 1), 
the cross section acquires the correction of the order k\ (we proceed as in Sec. 21. 4 
in [3] but taking into account polarization): 
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Here A c is the electron Compton wavelength, k 2 is the mean value of k 2 (q) for 
the atomic density n a (g) = exp(— g 2 / 2uf) / 2nul, the function h({u\/a 2 ) reflects a 
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nongomogeneity of atomic distribution in crystal, the function f(Za) represents 
the Coulomb corrections: 



r r e ~x/vi 2u 2 

K 2 = / K 2 (Q)n a (Q)d 2 Q = / K 2 (x)dx, 7]i = — 

J J m ai 

o s 



oo i 

/(0 = Re [V(l + - ^(1)] = e E 7 2 ■ . 2V (3-16) 

where V'(C) is the logarithmic derivative of the gamma function, Ei(— x) is the 
exponential integral function. It is seen from Eq. (|3.15|) that in a weak field limit 
(k\ <C 1) in the logarithmic approximation the integral positron polarization is 
C + = C C ^3/7. 

The probability of coherent pair creation becomes comparable with the Bethe- 
Maximon probability at K m ~ 1. In this situation Eq. (j3.12|) still has quite sat- 
isfactory accuracy (error < 15%), while Eq. (j3.15|) becomes inapplicable. The 
relative contribution of the LPM effect to the total probability of pair creation 
by a photon was found recently in [Hj. This contribution has the maximum 
about 5% at K m ~ 1 for heavy elements. In this paper we neglect the LPM 
effect. The differential cross section of incoherent pair creation in the effective 
crystalline field for polarized particles in the logarithmic approximation for an 
arbitrary photon energy can be presented in a form 

da p (£ + ) 2Z 2 a 3 °f ( x\ dx . ln(l + ^ 



/( x \ dx 
exp )f{x,y,€+) — , L(ki)=L u + 



dy 15m z J \ rji J r\\ 3 

(3.17) 

where 

fix, y, £+) = (1 + £+)fi(x, y) + [y 2 (l + £+) + (1 - yf{l - £+)] f 2 (x, y). (3.18) 
Here 

h(x, y) = z A T(z) - 3z 2 V(z) - z\ f 2 (x, y) = (z 4 + 3z)T(z) - 5z 2 V(z) - z\ 
z = z(x,y) = [k(x)u(1 - y)r 2/3 , y = — , (3.19) 

LO 

where the function k(x) is defined in Eq. fl3.8|) . Here T(z) is the Hardy function 

00 / 3\ 

T(z) = J sin [ ZT + Y)- ^ 3 - 2 °) 
o ^ ' 

When the influence of axis field on the incoherent process is weak [k\ < 1,2 > 
1), one can use decomposition (see Appendix D, Eq. (|D.3|0 

1 2 4D 64 86 

T(z)^- + - + ^, / 1 ( Z )~ 5 + H JJ / a ( z )~io+~ (3.21) 

z z 4 z' z 6 z 6 
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Substituting this decomposition in Eq. (|3.17| ) and integrating over y we obtain 
within the logarithmic accuracy Eq. (j3.15|) . 

Integrating Eq. (j3.17J) over the positron energy we find for the cross sections 
contained in Eq. (j3.14|) 



dx 



AlZ^ of* f f ( X \ r 

Gp = i5m 2 Lv<yKir ^ J dy J exp ( — ) lAfof) + 2 v 2 M x >y) 

171 oo V 

0£+ = 15m 2 L u( K rn) J dy j exp y-—J h( x iV)— ( 3 - 22 ) 

Here we introduced the functions L p (n m ) and Lg + (K m ) which refined the loga- 
rithmic approximation Eq. ()3.17j) for Ki < 1 and for K\ 3> 1- These functions are 
obtained by means of interpolation procedure (|14j): 

11 6 — + 3^ 
L p ( Km j - L u - - + - Jn , 



L ?+ ( Km ) - L u - - + - In g (3.23) 



In Fig. 2 the integral probability of pair creation by a photon in tungsten, axis 
< 111 >, T = 100 K as function of photon energy is shown. The curve 1 is the 
incoherent and the curve 2 is the coherent contribution while the curve 3 is their 
sum giving the total probability. In low energy region the LPM effect diminishes 
slightly the total probability (see ^1]). In this region the incoherent contribution 
dominates, these contributions are equal at uo ~ 10 GeV. At higher energies the 
coherent contribution dominates achieving maximum at to ~ 1.2 TeV, while the 
incoherent contribution aims for zero. 

In Fig. 3 the positron polarization vs relative positron energy e/uo is given: the 
curve 1 for uj = 12 GeV, the curve 2 for u = 22 GeV, the curve 3 for u = 100 GeV. 
Both coherent and incoherent contributions are taken into account. It is seen 
that shown dependence is approximately universal. When e/oj — > the positron 
polarization is tending to the Coulomb limit Q c = —1/3 since in this limit the 
coherent contribution is negligible. In the limit e/oj — > 1 the polarization £ + — » 1 
because of helicity transfer, which is common for both coherent and incoherent 
contributions. 

In Fig. 4 the integral positron polarization £ + as a function of photon energy 
uo is shown. Both coherent and incoherent contributions are taken into account. 
At uj — > this is is the Coulomb limit ( c = 3(1 — l/7L u )/7 (L u is defined 
in Eq. (|3.15j) ). At intermediate energies the shown dependence is the result of 
interplay of the coherent and the incoherent contributions, while in the high 
energy region we see the behavior of the coherent contribution. 

Now we turn to the situation when the photon energy is very high and influ- 
ence of crystalline field becomes very strong (ki ^> 1). In this case we have (see 
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Appendix B) for Ea. (f3~T4]) (a p is given in Eq. (21.34) in g|) 

8Z 2 « 3 r 3 (l/3) 9/ o , . , — ^ . , Z 2 o?LU x ) 
a p ~ 25m2 1 7 J 3 2 / 3 r(2/3)^-^L p (^) ~ 3.23 

*e + * ^ K - 4 /3 L Km ~ 3.95 — 3.24 

where the average ~k? is defined in Eq. (J3.16|) . It should be noted that the asymp- 
totic expansion of at K\ ^> 1 becomes valid only for very large values of K\ 
when the ratio a^ + ja v is very small. This is the consequence of specific form 
of the next term of decomposition of = + er|^ + . . . . The term aj^ 
(Eq. (|3.24j) ) is calculated in Appendix C. The next term has a form 

{2) Z 2 a M(lnKi + £) 
a l = —2 ~2 L (+ (n m ). (3.25) 



For example, for tungsten crystal, axis < 111 >, T= 100 K, the obtained in 
numerical calculation values of the coefficients are A ~ 10, B ~ —1.7. 

In the region k s 3> 1 (the constant field approximation is valid if k 1 J z '&q/'& v <C 
1) the integral probability of coherent pair creation (the terms oc £ + are calculated 
in Appendix B , the integral probability for the unpolarized case is given in 
Eq.(12.16) of jlj) is 

W?{£ + ) = c aV \ /3 [lnK s + B 1 (r ] ) 



+ ^+ f 1 ~ + c 2) J 1> (3-26) 



20 s+ V «2/ 8 

where (see Appendix B, Eq. (jB.5|0 

5 3^(2/3) Ci = ^Lim_ . 527i 

7vr 2 1 /sr(l/3) ' 3 19 / 6 T 3 (2/3) 

C2 = l n2 -^-C+^ 1.067; /3 - 7 



2 ' ' 1 + r/ 

BM = -0.374 - 3.075/3 2 / 3 (l + 1/9 + ^/? 2 ) + (2 + 5/J + g/3 2 ) (3.27) 

When k s = UiKx/a s 3> 1 the contribution of incoherent process is very small and 
can be neglected (see Eq.(21.35) of [I]): 
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4 Modified theory of coherent pair production 



The estimates of double sum in the phase A p \ made at the beginning of previous 
section: ~ (i9 v /$ ) 2 ty remain valid also for t? > $v, except that now the factor 
in the double sum is (^/^o) 2 < 1, so that the values |giir| > 1 contribute. We 
consider first the limiting case d v , then this factor is small and exp(iA p i) 

can be expanded accordingly. As a result Eq. (j2.9j) acquires the form 



„, rcoh , \ iam 2 de + 7 dr ( .mVr\ ^ G(q)G(q') /N 

dwg_\u) = — — — / — — ex P i — - 2 / (q±q±> 



m 2 uJT 



(Pp 2 (t,+) sin(g||r) sin^r) + i<p p i(£+) V(q\\, q[\ 

d 3 r 



exp[-z(q + q')r]. (4.1) 

The integration over coordinate r in Eq. (j4.1)l is elementary and gives 5q+q',0) 
after which the sum over q' and the integrals over r are easily calculated by 
means of the theory of residues. Finally we obtain 



^>) = ^£|G(q)| 2 ^ 



, . , . 2m 2 u ( ', , m 2 ui 



£_ 



■ £+q \\ 



m 2 uj 



2ss> 




For unpolarized photons (£ + = 0) Eq. (j4.2|) coincides with the result of standard 
theory of coherent pair production (CPP), see e.g. [TO] . 

In the case k s ^> 1 (k s is defined in Eq. f)3.8jl ). one can obtain from gen- 
eral expression Eq. (|2.1()j) the expression for spectral distribution, the region of 
applicability of which is broader than that of standard CPP theory. For this 
purpose it is necessary to take into account that the phase A p \ Eq. fl2.10|) has for 

1\\ + 9|i 7^ terms of the order (i9„/i?o) 3 / and (#u/$o) 4 / K l 2 ( K ' s — k, s e^e + /lu 2 ) 
which can be small even for $ < i9 v if k s ^> 1. Therefore, assuming that these 
contributions are small, we carry out the corresponding expansion of exp(iA p i), 
while the term with q\\ + q'» = in the double sum in A p i will be retained in the 
exponent. As a result we obtain an expression which coincides in the form with 
Eq. fl4.1|) where we must make the substitution 

exp \i— j - exp ^— j, m. = m (l + - ) . (4.3) 

Above the parameter g (Eq. (|2.1jl ) has the form 

ff 1 \^ r^^^S^s! [a a a 1 l g (q)l 2 qi /„ ^ 

/ m q,q' q W q \\ q,9||^0 m ?|| 
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sin(<7|| + q\y)r 

and in the term — ; — ^ — it is necessary to assume that on + q<< ^ 0. The 

(o||+gf|)r y " y " r 

remaining calculations are carried out in the same way as in the transition from 

Eq. ()4.1|) to Eq. ()4.2|) . The final result can be presented in a form 



Suj q g|| 



h g(y±jt) lfi T 

(2 + g)y r V r 



#(t -t), (4.5) 



where 



In derivation of Eq. (J4.5j) the higher order terms over Tg/r (2 + g) are omitted. 
Equation (J4.5|) is not more complicated than Eq. (|4.2j) but has a significantly 
broader range of applicability. 

The spectral distributions Eqs. (j4.2|) and ()4.5|) can be much higher than the 
Bethe-Heitler pair production distribution for small angles of incidence $o with 
respect to selected axis. For the case $0 <C 1 the quantity qn can be represented 
as 

The main contribution to the sum in Eqs. (j4.2j) and ()4.5)1 for small do is given q 
with m = 0, then 

/ 2tt 2tx \ 

q\\ ~ I — k cos ip + — Zsin<^J # , (4.8) 

where / and h are the characteristic periods of the potential in the plane trans- 
verse to the considered axis, cp is the angle of the projection n onto this plane 
with respect to one of the planes containing the selected axis, k and I are integers. 

Let us consider the spectral distribution of pair production in the extreme 
limit when the parameter s = 2uj\q\\\ m i n /m 2 ~ ujdQ/m 2 a s ^> 1. In this case the 
maximum of distribution is attained at such values of $0 where the standard 
CPP becomes inapplicable. Bearing in mind that if s 3> 1 and $0 ~ Vo/m than 
k s ~ s ^> 1, we can conveniently use the modified theory of CPP. Utilization 
of the modified theory for these values of s and g gives the exact position of 
maximum of the spectral distribution (in region where r ~ ro) and the value of 
the total probability within logarithmic accuracy (Ins ^> 1). 

The transverse component of the vector q in Eq. (j4.7|) can be selected in a 
such way that the spectral distribution given by Eq. (j4.5|) has a sharp maximum 
near the end of the spectrum at y m (l — y m ) = (2 + g)/2s with relatively narrow 
(in terms 1/s) width Ay ~(l + ^/2)/s 

V d y m.ar. 4 ( 2 + Q) V T J 2 + Q 
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It is seen that the maximum of the spectral distribution with the opposite helicity 
(£+ = —1) is suppressed as 1/tq. At r > r one have to take into account the 
next modes of q\\ (see Eq. (j4.7|0 . In this part of spectrum the suppression of the 
probability with opposite helicity is more strong, so the created positrons have 
nearly complete longitudinal polarization. 

Bearing in mind that Ay ~ (2 + g)/2s we find that for £ + = 1 

dN + ( dWj. + \ agm 2 g e e 

dt \ dy J max 8lj L rad u 

where L ra d is the radiation length in a corresponding amorphous medium, e e = 
m(lQirZ 2 a 2 n a \ 3 c L )-\ L = \n(183Z-^ 3 - f(Za)), the function f(Za)) is de- 
fined in Eq. ()3.16|) . N + is the number of created positrons, for tungsten one has 
e e ~ 2.5 TeV. Thus, the above analysis shows that the considered mechanism of 
creation of longitudinally polarized positrons is especially effective because there 
is a gain both in the monochromaticity and the total yield of polarized positrons 
with the energy e + ~ u>. 



5 Conclusion 

It is shown above that at high energy in the process of production of electron- 
positron pair with longitudinally polarized particles by the circularly polarized 
photon in an oriented crystal the phenomenon of helicity transfer takes place 
in the case when the final particle takes away nearly all energy of the photon. 
This is true in the constant field limit t9 ^ Vo/ m as we h as i n the coherent pair 
production region $ > Vo/m. 

In crossing channel: the radiation from longitudinally polarized high energy 
electrons in oriented crystals is circularly polarized (^ 2 ' — > 1) near the end of 
spectrum [TJ] also in both regions. This is once more the particular case of 
helicity transfer. 

So, the oriented crystal is a very effective device for helicity transfer from a 
photon to electron or positron and back from an electron to photon. Near the 
end of spectrum this is nearly 100% effect. 
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A Appendix 



Integral polarization of positron in constant field, asymp- 
totic 

It is instructive that from the expression for Ro(X) (Eqs. ()3.7|) . (|3.6jl ) one 
can find integral degree of positron longitudinal polarization in an external field 
characterized by parameter k for any u: this is ratio of coefficient at £+ and the 
rest part which enters in the probability for unpolarized particles 



Fx 



C = it, Fx = - 



F 

i 



dy 



y 



dy 



i-y 



UK 2/3 (X)+ J K 1/3 (z)dz 



i-y 



y 



K 2/3 (X) + / K 1/3 (z)dz 



3«y(l -y)' 
(A.l) 



where y = s + /uj. This result follows also from Eq.(3.70) [4 . The integrand in 
the first term of Fx is antisymmetric at substitution y <r-> (1 — y) and because of 
this doesn't contributes into integral polarization. 

Let us find the asymptotic values of C- At k <C 1 it is convenient to write the 
function F\ in the form 



F x 



/v / Kl/s(z)dz = I j7=rJ K ^ dz 



dx 



V a/(4y(X-y)) 



i \Jx(x- 1), 



(A.2) 



where a = 8/3k, x = l/(4y(l — y)). In the limit a ^> 1 one can use the standard 
expansion of Ky 3 (z) for z ^> 1 and the last integral in Eq. (jA.2|) becomes 



J K 1/3 (z)dz c ^ J e —dz 

ax ax 

Taking the integral over x we find 



Ti e~ 



2 wax 



Fx 



V2c 



while the function F Eq. (jA.l|) in this terms (see Eq.(3.58) in is 

3tt 



So, we find 



2y/2c 



(A.3) 



(A.4) 



(A.5) 



(A.6) 
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In the limit k> 1 we present the integral Fi as 



= / — / K 1/3 (z)dz = + F[ 2 \ a < 1, y < 1, - < 1 



* o/Jz/Cl-,/) 
1 , oo 



1 , oo 



3/0 0/42/(1-2/) 



7T 

73' 



2/0 , oo 2/0 , oo oo oo 

f[ 2) = J^ J K 1/3 {z)dz * J ^ J K 1/3 (z)dz= J ^-jK 1/3 (z)dz 



^ o/4»(l-i/) 



V 

3 a/42/ 



a/42/o 



7T 

71 



+ V 



+ ln2 



(A.7) 



where ip(x) is the logarithmic derivative of the gamma function. In calculation 
of F} 2) we substitute the variable s = a/Ay and than performed integration by 
parts. Finally we have for F\ in the limit k>1 



Using the two first terms of decomposition of total probability of pair creation 
(see footnote at p. 86 in 4 j) we obtain for integral probability of polarized pair 
creation: 



am 



2uo 



where 



5r(5/6)(2/3) 1/3 
D = — K /.J,l J — = 0.37961. 



(A.10) 



14T(7/6) 

So we have for integral degree of longitudinal polarization of created positrons 

In k- 1.126 



C(«) 



1.139k 2 / 3 -2" 



(A.11) 



B Appendix 

Coherent integral polarization term in a crystal at k s ^> 1 

Here we consider situation when the value of parameter k(x) on the boundary 
of cell is small: k(xq) ~ 2k s /x^ 2 <C 1 and one can extend the integration interval 
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over x up to infinity. The term oc £ + in the integral probability contains the 
integral (see Appendix A, Eq. (jA.2|) ) 



oo 1 , oo 







F ax = J dx J — J K 1/3 (z)dz, a(x) 



8 4(x+l)v^ 



3k,(x) 



(B.l) 



where g = a(x)/(4y(l — y)) The integration interval over x we split into two 
parts 



1. oc Oljqj 2* oc oc^^ 



where h?J 2 ^> x c ^> 1. In the first interval we have 

oo 1 oo 



Xc 



dx — dx 



to 



91 (as) 
1 oo 

dy 



1 , oo 

2 rtfy 



— / Ki/ 3 (z)dz — ~ — / K 1 / 3 (z)dz 



V 



3J y 





F (ll) , ^(12) , ^(13) 
ax ' ax ' i 



(B.2) 



where = x 3//2 /(3k s ?/(1 — ?/)). In the first integral over x : i* 1 ^ 1 ) = J dx . . . 

o 

we reverse the integration order, introduce the variable s = and integrate 

over s by parts. We get 

= (3.,) 2 / 3 / [y(l - y)f^ J s ^K 1/3 (s)ds = ^ (B.3) 



/ 3 ) 

In calculation of integrals over ?/ in the second interval over x, as well as in 
remaining integrals in Eq. (jB.2|) we use the results of previous Appendix. We 
have 



F (2) , F (12) , F (13) 

ax ' ax ' nrr 



7T 



In 



1 + x 



c&c H — In 



2 f\l Q 



c 



(B.4) 



Adding Ea. (|ET3j) and En.(jET4|) we find 



2/3 



2/3 



(In k 8 + c 2 ) 



g5/3 £3(2/3) 

2 4 / 3 r(l/3) 



2.295280.. 



Cl = 4?r2l 1 / l r iyf^ ^ 0.526820.., c 2 = ln2 - — - C + - ~ 1.066625.. (B.5) 
27- 3V6r 3 (2/3) 2 2 v ; 
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C Appendix 

Incoherent integral polarization term in crystal at Ki ^> 1 

At integration over y the term containing /^(a?, y) in Eq. ()3.18|) vanishes due to 
symmetry of integrand at y «-> 1 — y so the integral polarization term (containing 



oo 1 

e"^ 1 - / h{x, y)dy, 



(C.l) 



o 



Let us consider the integral 

1 1/2 



Fi = J h(x,y)dy = 2 J h(z{x,y))dy, z = [k V {\ - y)}- 2/ \ 



(C.2) 



where fi(z) is defined in Eq. (|3.19[) . Using the differential equation for Hardy 
function Eq. (jD.6|) we can represent fi(z) in the form 



f 1 (z) = z 3 T>(z)-3z 2 V(z). 
Let us split the integration interval in Eq. (jC.2|) into two parts 

1- < y < y , 2.y <y< 1/2; k' 1 < y < 1. 
In the first interval z ~ (/q/)~ 2 / 3 and the corresponding integral is 

dz 



(C.3) 





--J 


K J 

oo 


3 






Jdz 




.0 



| (z 3 T"(z) - 3z 2 V( 



,5/2 



V?T"(z) - ~-j=T(z) 



z 



(Ky )- 2 / 3 « 1(C4) 



Now we turn to the first integral 



OO / 

/ v zV ' [ 

n \ 



z) -j=T'(z) \ dz 



Using the representation Eq. (j3.2()j) we get 



oo . oo 



7 °° 

2/ T '< 2 



OO oo 



dz 



(C.5) 



T'(z) = [ —= [ t cos ( zt H ) dr = 2 / dx [ r cos ( x 2 t H \ dr 

K ' Jl J Jli I 3/ J J 1 3 I 

V V x 7 00 x 7 



2 J dx J v^cos ( x 2 + — j dr = —= J dy J dxcos(x 2 + y 



oo oo 





oo 







— = / rcos(r 2 W = — = / cossrfs = 0. 
y/3J \/3J 



y/3- 



(C.6) 
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The second integral in Eq. (jC.4j) is 

z / n \ z o 



= 2 -zT - 64 /2 T'(0). (C.7) 



In the second interval yo < y < 1/2 the variable z < 1 than the integral in 
Eq. dUl is 



if > ~ -2 

.W 



1/2 r 

3T'(0) 



K?y*(l-yf K 4 / 3 y 4 ^(l - V ) 4 / 3 



dy. (C.8) 



In the first term (oc k 2 ) we retain only the main term oal/yo and in the second 
term we carry out the integration by parts. We find 

p (2) 2_ _ 18T'(0) 18 ■ 2^T(0) _ 24T'(0) f dy 

Substituting Eq. ()C.5|) -Eq. ()C.7|) into Eq. (jC.4|) and combining the result with 
Ea. (jC.9j) we obtain 

Fl = F m + ff) ^ 2 J^1 ( 3 . 2 -V3 _ 6l ) , 

^l y^ty)^ *™"' F ^^- (C ' 10) 

Integrating in Eq. fjC.lj) over s = xjr)\ {k(s) ~ K\\/2s/{l + s)) and using the 
numerical value of the integral 

OO 

2- 2 / 3 j e- s s- 2 '\l + s)^ = 2.561.., (C.ll) 



we find 

D Appendix 

The Hardy function 

The Hardy function 



a u ~ 3.95 (C.12) 



2 4 /3 ' 



T(s) = J sin I + — I dr (D.l) 
o ^ ' 
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is encountered in the theory of electromagnetic processes in an external field. 
At z <C 1 the decomposition of T(z) is 

i f^ (-3^) fc r (k + i\ cos (k + i \ _ rW3) r(2/3) ^ 

1 W- 32/3 Z. fc , ^ 3 J C0S V 3 J " 2-3 2 /3 + 2-3V3 2; 2 + -' 

(D.2) 

and T(0) = 0.643950.., T'(0) = 0.469447... 

At z 3> 1 there is asymptotic series of T(z) over l/z 3 : 

_ N 1 ~ (3Jfe)! 1 1 / 2 40 2240 \ _ . 

For calculation it is convenient to use the following representation of the 
Hardy function and its derivative 



00 , 
)=/sin( 



\/3 n\ ( zt t 3 



T(z) = I sin ( — + - I exp I — — - — I dr, 



00 , 

V{z) = J cos l^zr + I ) exp f - £ ] rrfr. (D.4) 



V3 ' - - ! 





It can be obtained using the expression 



2 3 



T(z) = Im Jexp (^zt + y ^ dr (D.5) 



after a turn of integration line by the angle n/6. 
The Hardy function satisfies the equation 



T"(z) - zT(z) = -1. (D.6) 
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Figure captions 



• Fig.l The spectral probability of pair creation dw^ + /dy, the curves 1 and 
2 are for energy u;=22 GeV, the curves 3 and 4 are for energy c<j=100 GeV, 
the curves 5 and 6 are for energy cj=250 GeV. The curves 1, 3 and 5 are 
for £ = 1, and the curves 2, 4 and 6 are for £ = —1. 

• Fig. 2 The integral probability of pair creation by a photon in tungsten, 
axis < 111 >, T = 100 K as function of photon energy. The curve 1 is 
incoherent contribution, the curve 2 is coherent contribution the curve 3 is 
their sum giving the total probability. 

• Fig. 3 The positron polarization vs relative positron energy s/uj in tungsten, 
axis < 111 >,T = 100 K. The curve 1 for uo = 12 GeV, the curve 2 for 
uj = 22 GeV, the curve 3 for u = 100 GeV. Both coherent and incoherent 
contributions are taken into account. 

• Fig. 4. The integral positron polarization £ + as a function of photon energy 
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